DIAS-STP-10-01 



Chiral fermions and torsion in the early Universe 

Brian P. Dolan 
Dept. of Mathematical Physics, NUI Maynooth, Ireland 
and 

School of Theoretical Physics, 
Dublin Institute for Advanced Studies, 
10 Burlington Rd., Dublin 8, Ireland* 



o 

(N 

< 

(N 

u 

cr, 



> 
m 



On 
O 



X 



Torsion arising from fermionic matter in the Einstein-Cartan formulation of general relativity is 
considered in the context of Robertson-Walker geometries and the early Universe. An ambiguity 
in the way torsion arising from hot fermionic matter in chiral models should be implemented is 
highlighted and discussed. In one interpretation, chemical potentials in chiral models can contribute 
to the Friedmann equation and give a negative contribution to the energy density. 
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It is often the case that quantum matter acts as a 
source for a classical field in situations where quantum 
aspects of the field itself can be ignored. This approxima- 
tion has proven extremely useful for Einstein's equations 
where 



Gab — 8ttG < T ab >, 



(1) 



works well when the matter source is degenerate 
fermionic matter, where <> is a quantum expectation 
value, and for thermal radiation, where <> is a ther- 
mal average of photons. There are difficulties with this 
approach however, not least that the singularities inher- 
ent in fully fledged quantum field theory for the sources 
render (1) ambiguous and some criterion for cutting off 
the integrals must be introduced. For example it is well 
known that a naive calculation of the vacuum energy den- 
sity of the standard model of particle physics leads to far 
too high a value of the cosmological constant to be com- 
patible with observations [1]. Nevertheless (1) seems to 
work well in the early Universe when the dynamics is 
dominated by radiation, as long as temperatures are well 
below the Planck temperature. In the radiation dom- 
inated Universe Einstein's equations boil down to the 
Friedmann-Robertson- Walker (FRW) equation, ignoring 
spatial curvature this is 



a 2 8nG m Ar 4tt 3 T 4 

- 2 = — < T 00 >= N eff 



45 m 2 pl 



(2) 



where a(t) is the cosmological scale factor and N e f / is the 
effective number of degrees of freedom in the relativistic 

gas, 



N eff :=N B + 



N-) 



(3) 



with Nb the number of bosonic degrees of freedom (2 
for photons), and N + and iV_ the number of positive 



chirality and negative chirality fermionic degrees of free- 
dom respectively (for a standard model neutrino N + = 2, 

The 



N. 



0: for a Dirac fermion 7V 4 



2), [2]. 



Planck mass, m Pl = G (we use units with % = c = 1), 
appears in (2) not because we are considering a theory 
of quantum gravity but because of the quantum nature 
of the source for classical gravity. 

In the Einstein-Cartan formulation of general relativity 
fermionic matter is expected to induce torsion (recent 
bounds on the magnitude of torsion have been derived 
from tests of violation of Lorentz invariance [3] and from 
cosmic microwave polarization [4]). When the connection 
is varied in the Einstein-Cartan action the torsion two- 
forms T a = ^T a oc e b A e c are determined by a spinor field 

via the algebraic equation 



r a = 2nGe abcd {y 1 5 1 d y)e b A e c , 



(4) 



a, b, c, . . . are orthonormal indices (for a review of torsion 
in Einstein-Cartan formulation in general see [5] and [6]). 

In the spirit of (1) the equation of motion (4) would 
be interpreted as 

Ta.bc = 4nGe abcd < vFyV* >= -AirGe abcd < jf >, 

(5) 

where jg = ^>^ a ^ 5 ^ is the axial current. As is well known 
fermions generate torsion in the anti-symmetric class of 
tensors, according to the classification of [7]. 

We shall examine the effect of torsion arising from rela- 
tivistic fermions in the early Universe, assuming isotropy 
and spatial homogeneity of both the geometry and the 
matter. It will be assumed that the metric of Robertson- 
Walker type and that the energy-momentum is of the 
form 



Tab = 



P 
PS, 



(6) 



where the density p and pressure P are homogeneous 



2 



and depend only on time and i,j — 1,2,3 are space- like 
indices. 

The Riemann tensor involves the square of the connec- 
tion and the net effect of including the torsion (4) into 
the gravitational connection is that Einstein's equations 
are modified to 



8irGp 



2a a~ t" 

+ — = 8ttGP, 



(7) 
(8) 



i 2 - 

a a 2 ' 4 

where a(t) is the Robertson- Walker scale factor and 

r 2 = -T a r a = -167r 2 G 2 (*7 5 7 a *)(*7 5 7 a *) (9) 

(the metric signature is (— , +, +, +), the Clifford algebra 
convention is {7 a ,7 6 } = —2r] ab , with 7 hermitian, and 
spatial curavature is taken to be zero). Equations (7) 
and (8) are not independent and are related by the (first) 
Bianchi identity which is analysed below. Eliminating r 
gives the usual relation between the acceleration and the 
density and pressure, 



a 
a 



4ttG 



(p + 3P). 



(10) 



Since fermions constitute quantum matter, it seems 
natural to interpret (4) in the early Universe as meaning 
a thermal average (5). However there is an ambiguity as 
to whether (9) should be interpreted using 



< (*7V*)(*7 7«*) >=< j 5 a j 5 , a > 



(11) 



< * 7 5 7 a * >< *7 5 7a* >=< jS >< jh.a > (12) 

These are different in general. The former can be Fierz 
re-arranged to give 



< (*7 5 7 a *)(*7 5 7 a *) > 



4 < (* 1 f *_)(* t _*+) > 



(13) 



where 1 $>+ and are the positive and negative chiral- 
ity components of \Jj. In a first quantised theory this is 
positive definite for Dirac spinors and vanishes for Weyl 
spinors [5], hence r 2 < in (9), making the torsion space- 
like (in the classification of [7] this is denoted As). The 
cosmological consequences of this formulation in infla- 
tionary models are explored in [8] . The same philosophy, 
applied to spin densities rather than the pseudo-vector 
vj/^s^a,^ is f n owcc i i n [9] anc i [iq]. 

We reach a radically different conclusion if we use (12), 
which follows from taking the thermal average of (4) be- 
fore calculating the Riemann tensor. Applying the usual 
Robertson- Walker assumptions of spatial homogeneity 
and isotropy to the connection, and hence the torsion, 
we would conclude that, in the cosmic frame, 



< 3l >= 0, 



(14) 



while 



< is >= ( n + - n +) - ( n - - n -) 



(15) 



where n+, n+, n_ and n_ arc thermal averages of the 
number densities of positive chirality particles, positive 
chirality anti-particles, negative chirality particles and 
negative chirality anti-particles respectively. In contrast 
the vector current j a = \D r 7 \E f ) has 



< j° >= (n + - n+) + (n_ - n_) = n — n, 



(16) 



where n = n + + n_ is the number density of all particles 



and n 



n- the number density of anti-particles. 



The integral of j° over a co-moving 3-dimensional volume 
V gives the total number of fermions in that volume, 
counting anti-particles as negative, 



j°a 3 (t)d 3 x = N-7J, 



(17) 



where J\f = J v na (t)d x, etc. The corresponding num- 
ber for the axial current, 



X 



j%a 3 (t)d 3 x = 7V+ - Af- - (AT+ - JJ-) := N A (18) 



we shall call the axial particle number in the volume V . 

In the abscence of chemical potentials, for a thermal 
distribution with the temperature much greater than any 
particle masses, 



n± = n± = — 



3C(3) 



4 ir 2 



rp6 



(19) 



so JV = N and the torsion vanishes. A non-zero chem- 
ical potential is necessary for any asymmetry between 
particle and anti-particle numbers. We see here that a 
chemical potential can also generate torsion, since then 



n± — n± 



iu ± T 2 + of^V 



G 



\ T J 



(20) 



for fi± « T, where p + and are chemical potentials 
for positive and negative chirality particles, see e.g. [2]. 
In a chiral theory, such as the Standard Model, p + and 
/U_ can be different in general and 



<j!>=-(p+-v-)t 2 , 



(21) 



where we have ignored terms o (7^) • In a 
Robertson- Walker Universe undergoing adiabatic expan- 
sion, a oc 1/T, the thermal average of the axial current 
V a < j'5 >= is conserved, and hence Ma is constant, 
if and only if (/i + — /i_) oc T is linear in T. The to- 
tal fermion number, JV — IV, is constant if and only if 
+ H-) oc T (particle masses are being ignored here). 
To summarise, in general equations (5) and (12) give 



T 2 = 



4^ 2 G 2 



(22) 



3 



which is positive in any chiral model for matter with 
/!+ 7^ For a model with Nf different types of 

fcrmion each fermions species can have different chem- 

(k) 

ical potentials p^ , where we label the species with an 
integer k, and then p± is always understood below to 
mean p± = J2k=i f 1 ^' '■ Since positive and negative chi- 
rality particles in a chiral model can have different weak 
charges one expects that ^ p^f 1 in general, and the 
torsion can be non-zero, at least for the era before the 
electro- weak phase transition, [11]. 
Define 



r=—(p+-p-)T 2 , 



(23) 



in terms of which the non- vanishing components of the 
torsion are 



(24) 



(in the classification of [7] this is time- like, At). Rota- 
tional invariance of the thermal average is not incompat- 
ible with the conclusion of [12], where classical solutions 
of the Weyl equation were analyzed in spherical sym- 
metric space-times with torsion — thermal averages do 
not necessarily have the same symmetries as solutions of 
the equations of motion. The general form of the torsion 
compatible with Robertson- Walker symmetries was given 
in [13]. The fact that chiral fermions can have interest- 
ing consequences when torsion is taken into consideration 
was noted in the context of anomalies for lepton currents 
in the Standard Model of particle physics in [14]. 

Both (11) and (12) have interesting, though very differ- 
ent, cosmological consequences. The form (12), being the 
square of a vector, has a dual description as the square of 
a 3-form and as such is in the class of models described in 
[15]. Indeed a term of this form is present in the Landau- 
Ginsparg models discussed in [15], though the stabilising 
quartic term is absent and there is no kinetic term here. 
A kinetic term would require time derivatives of the tor- 
sion and so would go beyond Einstein-Cartan theory — 
such terms would be expected to appear in an effective 
action description of gravity involving higher derivatives 
and powers of the Riemann tensor. 

So which should one use (11) or (12)? Weinberg [16] 
takes the point of view that there is nothing special about 
torsion: it is just another tensor and one can always move 
it to the right hand side of Einstein's equations and con- 
sider it to be part of the matter rather than part of the 
geometry. We see here that, in the context of (1), there 
is an ambiguity. If the torsion terms are absorbed into 
the energy momentum tensor before expectation values 
are taken then it would seem that (11) is appropriate. 
In the Einstein-Cartan formulation however the torsion 
is determined by the equation of motion (4), in which 
the square of the torsion does not appear. If the gravi- 
tational field itself is not quantised, it is hard to see any 



interpretation of the equation of motion (4) other than 
(5). When the Riemann tensor is calculated it is then 
(12) that arises and not (11). Much of the literature has 
focused on (11), in a cosmological context for example 
(11) was used in [8]. In this paper the consequences of 
(5) and (12) will be explored and developed. 

We shall see that, in the context of the early Universe, 
the torsion can give a negative contribution to the en- 
ergy density. The mechanism here is different to torsion 
induced avoidance of the initial singularity due to spin 
fluids considered previously, [17-21], in which the spin 
density necessarily breaks either rotational or translation 
invariance. 

When there is torsion the Bianchi identity does not 
require that G ab be co-variantly constant, in general one 
has 



\7 b G ba = -T c bc G ba + -R adbc 



Td,bc, 



(25) 



where R adbc := je ada ' d ' R a , d , b , c ,e bcb ' c ' . In the case of 
Fricdmann-Robertson- Walker (FRW) Universes under 
study here only the second term on the right hand side 
contributes giving 



v ' gm - -iii<™>. 



bi 



0. 



(26) 



One strategy is to demand V b G ba = and use this to 
determine the torsion, implying that r « 1/a [22], but 
this is too restrictive for our purposes. Instead we take 
thermodynamic averages as above and use (23) for the 
form of the torsion. 

Assuming adiabatic expansion T oc 1/a, [2], a can be 
eliminated from Einstein's equations in favour of T to 



give 



f 

T 



8irG ir 2 
A-kG 



{bp + 3P) 



(,i + - M _) 2 T 4 

2tt 2 G 2 , 

-(P+ 



(27) 
2 T 4 . (28) 



3 v ' ' 9 

These two equations are not independent, the Bianchi 
identity (26) gives an equation relating p, P, T and p±. 
Expressing time derivatives as temperature derivatives, 
3j = T-gp, differentiating (27) and using (28) to eliminate 
T gives 

I^^-^V} -£(»-!■*), (29, 

where ft = p + P is the enthalpy density. For standard 
equations of state (relativistic gas, dust, cosmological 
constant) the right hand side of (29) vanishes, so the tor- 
sion described here necessarily requires a modification of 
the equation of state. Let p = p + Ap, P = P + AP and 
h = h + Aft, where po etc., satisfy a standard equation 
of state. Then 

^<<<"-> 2r2 } = §( 3a '-^)- < 30 » 



4 



For a relativistic gas with the full energy-momentum 
tensor traceless AP = 4p and (30) reduces to 

Two special cases are: 

• If axial particle number is conserved, — H-) = 
bT with b a constant, so Ap must have a T 6 com- 
ponent, 



4 _ nGb' 6 
12 ' 



(32) 



where the T 4 term is assumed to have the usual 
form for massless particles with N e ff degrees of 
freedom. In particular the torsion necessarily gives 
a negative contribution to the energy density in this 
case, though in realistic models this is a small ef- 
fect. We must take b « 1 to be consistent with 
our assumption that (/x+ — fi-) «T. This is com- 
patible with the observed value of the current ratio 
of Baryon to photon densities 



12C(3) 



10" 



(33) 



with riy = ^|^T 3 . Tj n i ess /U + and p,_ have oppo- 
site sign and the smallness of this number is due to 
a delicate cancellation between two large numbers, 
this requires b to be of the order of, or less than, 
10~ 9 at the present time. Then the ratio of the two 

terms in (32) is of order which is very 

small for T « mpi (obviously the analysis here is 
only valid for T « mpi where quantum gravity 
effects are assumed to be very small). 
This behaviour t ~ 1/a 3 for torsion arising from 
spin in the early Universe has been studied before, 
[20, 21]. 

• For any period during which p + — fj,_ is constant, 
independent of temperature, equation (30) gives a 
logarithmic correction to the the energy density, 



7r , irG . . o A ( T 

^30^ T4 -T2-^-^ Tln U 



(34) 



where T is an arbitrary constant that can be ab- 
sorbed into the definition of N e ff. Axial particle 
number is not conserved in this case and one ex- 
pects Baryon/Lepton number violating processes 
unless p,+ = — fi-. 

Of course the generation of Baryon or Lepton number 
requires P and CP violation as well as a period out of 
thermal equilibrium, according to the Sakharov condi- 
tions, [25]. The details would depend on the particular 



chiral model generating the chemical potentials and tak- 
ing the Universe through a period during which it is out 
of thermal equilibrium. 

At first sight if might seem disconcerting that energy- 
momentum does not appear to be conserved in this for- 
malism — because of (26) and the Einstein equations 
T a b cannot be co-variantly constant unless ar is constant. 
However an "improved" energy-momentum tensor, which 
is conserved, can be defined. We make the co- variant de- 
composition of the Einstein tensor 



Gab = Gab + AG ab 



(35) 



where Gab is the Einstein tensor constructed from the 
torsion-free connection. We similarly decompose the con- 
nection one-forms as 



uj a b =lj a b + Auj a b 



(36) 



with uj a b the torsion- free connection. Expanding Aw";, = 
A(jj a b, c e c the components Auj a b. c , being the difference 

of two connections, constitute a tensor field so (36) is 

o 

again a co- variant decomposition. Gab is the zero torsion 
Einstein tensor for which the first Bianchi identity implies 



o o , 

VbG ba = 



(37) 



o o 
where Vfc is the co- variant derivative using uj a b . From 

this follows 

V b G ha =V b (AG ba ) + Aoj b c . b G ca + Aoj a Ctb G bc . (38) 
We also have, by definition, 

V b T ba =Vb T ba + Aio b Ctb T ca + Acv a c , b T bc (39) 
for T ab . Einstein equations, G ab = 8irGT ab , now imply 

(40) 

An "improved" energy-momentum tensor can be defined 



Vfc (AG ba ) = 8ttG Vfc T ba . 



q-ab rpab 



8ttG 



AG 



ab 



(41) 



which is conserved using the torsion free connection, 
o 



VfcT' 



bet 



o. 



(42) 



For example if r = ( 27r 3 Gb )r 3 in a radiation domi- 
nated Universe, the improved energy-momentum tensor 
for FRW space-time with torsion is 



Tab = 



Po 

U 3 Oij 



b 2 irG 
24 



T 6 S ab 



(43) 



■5 



with po = ^N e ffT 4 . In fact both terms in (43) arc 
separately conserved with the torsion-free connection. 

Finally we observe that the geometrical significance 
of non-zero r follows from the anti-symmetrised action 
of two co- variant derivatives on an arbitrary vector field 
with components U a , 

[V a ,W b ]U c = -T d ab V d U c + R c dab U d . (44) 

In addition to the algebraic (rotation) term involving the 
Riemmann tensor there is a derivative term involving the 
torsion — a deficit displacement implying that parallelo- 
grams generated by parallel transport do not close. The 
deficit displacement in Robertson- Walker space-time de- 
scribed here is compatible with 3-dimensional rotational 
symmetry — a vector field with Robertson- Walker sym- 
metries must have U l = in which case 

[V l ,V J ]U° = -re^VkU (45) 

and this vanishes if U° independent of position. For any 
field compatible with the Robertson- Walker symmetries 
space-like parallelograms close with the torsion studied 
here. However they need not close for fields that do not 
share the symmetries of the background metric. 

The author wishes to thank the Perimeter Institute, 
Waterloo, Canada, and the Dublin Institute for Ad- 
vanced Studies for financial support. 



* Electronic address: bdolEm@thphys.nuim.ie 
[1] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1 . 
[2] S. Weinberg, Gravitation and Cosmology, Wiley (1972). 
[3] V.A. Kostelecky, N. Russell and J.D. Tasson, 
Phys. Rev. Lett. 100 (2008) 111102, arXiv:0712.4393[gr- 
qc]. 



[4] Moumita Das, Subhendra Mohanty and A.R. Prasanna, 
Constraints on background torsion from birefringence of 
CMB polarization, arXiv:0908.0629[astro-ph.CO]. 

[5] F.W. Hehl, P. von der Heyde and G.D. Kerlick, Rev. 
Mod. Phys. 48 (1976) 393. 

[6] I.L. Shapiro, Phys. Rep. 357 (2002) 113, arXiv:hep- 
th/0103093. 

[7] S. Capozziello, G. Lambiase and C. Stornaiolo 

Ann. Phys. 10 (2001) 713. 
[8] T. Watanabe, Dirac-field model of inflation in Einstein- 

Cartan theory, arXiv:0902.1392[astro-ph.CO]. 
[9] I.S. Nurgaliev and W.N. Ponomariev, Phys. Lett. 130B 

(1983) 378. 

[10] M. Gasperini, Phys. Rev. Lett. 56 (1986) 2873. 
[11] J.A. Harvey and M.S. Turner, Phys. Rev. D42 (1990) 
3344. 

[12] B. Kuchowicz, Phys. Lett. 50A (1974) 267. 
[13] M. Tsamparlis, Phys. Lett. 75A (1979) 27. 
[14] A. Dobado and A.L. Morato, Phys. Rev. D54 (1996) 
5185. 

[15] Inflation and dark energy from three-forms, T.S. Koivisto 
and N.J. Nunes, arXiv:0908.0920[astro-ph.CO]. 

[16] S. Weinberg, Physics Today, 59 April (2006) 15. 

[17] W. Kopczyfiski, Phys. Lett. 39A (1972) 219; ibid. 43A 
(1973) 63. 

[18] A. Trautman, Nature (Phys. Sci.), 242 (1973) 7. 

[19] G. de Berredo-Piexoto and E.A. de Freitas, On the tor- 
sion effects of a relativistic spin fluid in early cosmology, 
arXiv:0902.4025[gr-qc]. 

[20] M. Gasperini, Repulsive gravity in the very early Uni- 
verse, arXiv:gr-qc/9805060. 

[21] L.C. Garcia de Andrade and R.O. Ramos, Spin-Torsion 
in Chaotic Inflation, arXiv:gr-qc/9910053. 

[22] P. Minkowski, Phys. Lett. 173B (1986) 247. 

[23] M.C. Gonzalez- Garcia and M. Maltoni, Phys. Rep. 460 
(2008) 1. 

[24] C. Barragan, G.J. Olmo and H. Sanchis-Alepuz, 

Phys. Rev. D80 (2009) 024016. 
[25] A.D. Sakharov, JETP Lett. 5 (1967) 24. 



